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Scientific  Progress  and  Accomplishments  since  September  18  2013: 

In  this  final  report  we  will  present  the  major  aehievements  obtained  during  the  period  of 
duration  of  the  grant.  The  findings  here  presented  are  eurrently  being  summarized  in  a 
manuseript  that  will  be  eompleted  and  submitted  for  publieation  in  the  following  weeks. 


1 .  A  framework  to  study  eompetition  for  eellular  resourees  in  Escherichia  coli 

One  of  the  priorities  of  the  projeet  was  to  build  a  system  to  quantify  the  extent  of  the 
eompetition  of  different  genes  for  the  same  eellular  resourees.  We  generated  a  synthetie 
eireuit  (MBP-1.0)  that  allows  us  to  study  the  simultaneous  expression  of  two  reporter  genes 
eneoding  a  green  (GFP)  and  a  red  (RFP)  fluoreseent  protein.  The  eireuit  is  eontained  in  a 
plasmid  with  a  tunable  eopy  number.  One  of  the  reporters,  GFP,  is  expressed  eonstitutively 
whereas  the  other  is  under  the  eontrol  of  the  transeriptional  aetivator  LuxR.  When  the  input 
of  LuxR,  the  indueer  aeyl  homoserine  laetone  (AHL),  is  added  to  the  eulture  medium, 

LuxR  binds  to  the  promoter  Plux  and  RFP  is  produeed  using  resourees  that  were  initially 
alloeated  into  the  GFP  synthesis.  As  a  eonsequenee,  the  amount  of  GFP  generated  is 
deereased  (Fig.  1  A). 

We  monitored  in  a  flow  eytometer  the  produetion  of  both  reporters  during  exponential  eell 
growth  on  glueose  0.4%  as  the  sole  earbon  souree.  The  data  shows  that  as  the  eoneentration 
of  RFP  inereases,  the  eoneentration  of  GFP  deereases  by  up  to  a  69.4%  for  maximal  AHL 
(Fig.  IB,  left  panel).  This  observation  eorrelates  well  with  numerieal  simulations  (Fig.  IB, 
right  panel)  of  a  system  of  ordinary  differential  equations  that  aeeounts  for  the  eonservation 
of  ribosomes  and  RNA  polymerase  along  with  the  produetion  and  degradation  of  mRNAs 
and  proteins.  The  parameters  used  in  the  model  were  obtained  from  the  literature. 

We  eonsidered  several  alternative  reasons  to  eompetition  for  eellular  resourees  that  eould 
explain  our  observations.  All  of  them  were  disearded  experimentally  being  the  eontrols 
shown  in  Fig.  1C  the  most  representative.  Toxie  effeets  of  the  indueer  and  RFP  produetion 
were  ruled  out  using  a  eontrol  that  does  not  eontain  RFP  (MBP-dRFP)  and  a  eontrol  that 
produees  a  protein  different  from  RFP  (MBP-gapA),  respeetively.  GapA  eneodes  for  the 
glyeeraldehyde  dehydrogenase  from  E.  coli  and  it  is  very  abundant  in  the  baeterial 
eytoplasm.  The  eontrol  MBP-dRFP  does  not  display  an  effeet  on  GFP  whereas  the  eontrol 
MBP-gapA  shows  a  deerease  in  GFP  very  similar  to  the  original  in  the  eireuit  MBP-l.O 
(Fig.  ID).  As  a  eonelusion  of  the  experiment  we  demonstrate  that  it  is  not  the  espeeifie 
produetion  of  RFP  but  the  proeess  of  produetion  of  a  protein  what  is  deereasing  the 
synthesis  of  GFP,  due  to  the  eompetition  for  limited  eellular  resourees. 
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Fig.  1.  Schematic  representation  of  the  circuits  used  to  study  competition  (A).  GFP  is  constitutively 
expressed  and  RFP  is  under  the  control  of  activator  LuxR  and  input  AHL.  Curved  arrows, 
promoters;  hairpins,  terminators.  RFP  (upper  row)  and  GFP  (lower  row)  expression  dynamics  (B). 
Each  colored  line  corresponds  to  a  concentration  of  inducer  AHL  in  the  experimental  results  (left 
panel)  and  simulations  (right  panel).  (C)  Control  circuit  MBP-dRFP  does  not  contain  RFP  whereas 
in  MBP-gapA  the  RFP  gene  has  been  replaced  by  the  glyceraldehyde  dehydrogenase  encoding  gene 
(gapA)  from  E.  coli.  (D)  Dose  response  plots  for  circuit  MBP-1.0  and  the  controls  MBP-dRFP  and 
MBP-gapA  in  the  steady  state.  Left  panel  shows  GFP  and  right  panel  shows  RFP  (note  that  the 
controls  do  not  produce  RFP).  All  plots  represent  mean  values  and  standard  deviations  of 
populations  analyzed  by  flow  cytometry  in  three  independent  experiments. 


2.  Influence  of  RBS  strength  and  copy  number 

In  order  to  determine  the  relative  contribution  of  RNA  polymerase  depletion  and  ribosome 
depletion  to  the  GFP  decrease,  we  created  a  set  of  circuits  with  progressively  weaker  RBS 
strength  for  the  RFP  gene.  The  weakest  RBS  extents  minimal  demand  for  ribosomes  by 
RFP  mRNA,  providing  a  way  to  assess  the  contribution  of  depletion  RNA  polymerase 
depletion  on  the  reduction  of  GFP.  We  used  a  set  of  RBS  sequences  that  range  from  very 
strong  (MBP-1.0)  to  very  weak  undetectable  translation  of  RFP  (MBP-0.006).  The  dose 
response  curves  show  that  weaker  RBS  strengths  have  reduced  effects  on  GFP,  with  no 
appreciable  effect  in  the  case  of  MBP-0.006  (Fig.  2A).  Numerical  simulations  of  the  ODEs 
are  consistent  with  the  experimental  data  (Fig.  2B).  This  result  supports  the  hypothesis  that 
the  pool  of  RNA  polymerase  is  not  appreciably  reduced  due  to  transcription  of  the  RFP 
gene  and  hence  that  ribosome  depletion  is  the  contributor  to  the  68.4%  reduction  of  GFP 
observed  in  Fig  1. 

Another  way  to  reduce  the  coupling  between  RFP  and  GFP  expression  is  to  decrease  the 
plasmid  copy  number.  We  tested  two  lower  plasmid  copy  numbers  to  analyze  the  extent  to 
which  the  coupling  could  be  reduced.  Instead  of  replacing  the  origin  of  replication,  which 
may  in  turn  generate  artifacts  due  to  the  involvement  of  different  replication  machinery  for 
each  of  the  origins,  we  used  the  DIAL  strains  (Kittleson  et  al,  J.  Biol.  Eng.  201 1,  5:10).  Eor 
these  experiments  we  always  used  the  original  J92001  construct  but  changed  the  host,  to 
tune  the  number  of  copies  of  the  circuit.  The  nominal  values  in  the  exponential  phase  of  the 
copies  tested  were  64±2  (J,  high),  30±14  (H,  mid)  and  4±1  (E,  low).  The  results  show  that 
the  extent  of  competition  depends  on  the  number  of  copies  of  the  plasmid:  it  is  reduced  to 
55%  with  medium  plasmid  copy  number  and  to  a  29%  with  the  lowest  copy  number  tested 
(Eigs.  2B).  This  observation  is  in  good  agreement  with  the  model  prediction. 
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Fig.  2.  AHL  dose  response  plots  as  a  function  of  RFP’s  RBS  strength  and  plasmid  copy 
number.  Left  and  right  columns  show,  respectively,  experimental  and  computational  results  for 
both  RFP  (odd  rows)  and  GFP  (even  rows).  (A)  RBS  strength  effects.  The  numbers  following  the 
suffix  MBP-  indicate  the  relative  strength  of  the  RBS  for  RFP  compared  to  MBP-1.0.  (B)  Plasmid 
copy  number.  The  plasmid  MBP-1.0  was  tested  in  the  DIAL  hosts  JTK60  J  (61-65  copies),  H  (16- 
44  copies)  and  E  (3-5  copies).  These  copy  numbers  lead  to  68.4%,  28%,  and  29%  change  in  GFP, 
respectively.  Experimental  values  depict  mean  and  standard  deviation  of  the  populations  for  three 
independent  steady  state  experiments. 
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3.  Analytical  model  and  linear  relationships  between  reporters  coneentration 

To  diseover  how  key  parameters  eontrol  the  extent  of  eoupling  we  eonstrueted  an  analytieal 
model  that  prediets  how  the  coneentration  of  one  protein  (GFP,  for  example)  is  affeeted  by 
induction  of  a  different  protein  (RFP,  for  example)  due  to  sharing  RNA  polymerase  and 
ribosomes.  This  model  is  eonstrueted  starting  from  the  eomputational  model  used  for  the 
simulation  results  and  making  some  approximations  that  are  supported  by  the  data. 
Speeifically,  these  approximations  assume  that  the  eoneentration  of  free  RNA  polymerase 
and  of  free  ribosomes  is  suffieiently  smaller  compared  to  the  respeetive  dissoeiation 
eonstants  with  DNA  and  RBS,  respeetively.  The  opposite  limit  cases  when  some 
dissoeiation  eonstant  is  mueh  smaller  than  the  free  amount  of  RNA  polymerase  or 
ribosome  are  not  supported  by  data. 

Based  on  this  model,  the  eoneentration  of  protein  Pi  (i=l,2)  is  given  by 


AB  +  2d,[£,(A/^,)  +  (x,  /.fi:,)(m,  /!,)] 

i=i 

in  whieh  A  and  B  are  constants  independent  of  the  synthetie  eireuit,  xt  andyr  are  the 
available  eoneentrations  of  RNA  polymerase  and  ribosomes,  ATs  the  dissoeiation  constant 

of  RNA  polymerase  from  the  promoter  of  gene  i,  is  the  dissoeiation  eonstant  of 

ribosomes  from  the  RBS  of  gene  i  (inversely  proportional  to  the  RBS  strength),  Di  is  the 
DNA  eopy  number,  a,  is  a  eonstant,  m,  is  the  amount  of  mRNA  per  gene  eopy,  and  Si  is  a 
number  between  0  and  1  and  represents  the  level  of  induetion  of  the  promoter  (See  SI  for 
more  details).  This  expression  quantitatively  predicts  that  as  protein  1  is  indueed,  protein 
2’s  eoneentration  is  redueed  and  the  extent  of  reduetion  depends  on  the  various  parameters 
as  indieated  by  the  formula.  Speoifieally,  as  Ky  is  inereased  (weaker  RBS  strength  for  gene 
1),  the  reduetion  in  protein  2’s  eoneentration  is  less  severe.  Similarly,  as  the  DNA  amount 
is  deereased,  less  reduetion  of  protein  2’s  eoneentration  is  observed,  but  this  is  to  the 
expense  of  a  lower  value  of  the  eoneentration  of  protein  1 .  By  This  formula,  one  ean 
predict  that  a  linear  relationship  between  Pi  and  P2  holds; 


P2  =  -aPy  +  q 


in  whieh  a>0,  that  is,  the  eoneentration  of  the  two  proteins  is  eonstrained  on  a  linear 
manifold  when  protein  1  is  indueed.  The  slope  of  the  line  a  inereases  as  the  dissoeiation 
eonstant  Ky  is  deereased,  that  is,  the  RBS  for  gene  1  is  weakened.  The  intereept  q,  instead, 
stays  eonstant  as  the  RBS  strength  of  gene  1  is  ehanged  (Fig  1  A). 
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The  analytical  model  can  be  used  to  study  the  case  in  which  there  are  two  inducible 
modules  (Fig.  3B).  Specifically,  the  linear  relationship  for  protein  production  can  be 
applied  to  both  modules,  rendering  two  lines  with  different  slopes  and  offsets.  In  the 
absence  of  competition,  the  amount  of  proteins  produced  by  the  circuit  could  lay  anywhere 
in  the  surface  delimited  by  maximum  production  of  each  protein  (Fig.  3B,  left  panel). 
However,  because  of  the  competition  for  RNAP  and  ribosomes,  the  achievable  region  for 
protein  production  is  limited  to  the  area  comprised  by  the  linear  manifolds  of  both  modules 
and  the  axes  (Fig.  3B,  right  panel). 


mean  RFP  (AU) 


Fig.  3.  (A)  Linear  relationships  between  GFP  and  RFP  production.  This  plot  is  a  summary  of  the 
information  presented  in  the  rest  of  the  figures  and  compares  the  emission  of  GFP  and  RFP 
controlled  by  the  various  RBSs.  Experimental  results  (left  column  in  A  and  B)  are  compared  to 
computational  simulations  (right  column  in  A  and  B).  (B)  Regions  allowed  for  protein  production 
in  the  case  of  using  two  inducible  modules.  The  left  panel  represents  the  case  where  there  is  no 
competition  between  the  genes  and  each  of  the  proteins  can  be  produced  up  to  its  maximum  rate. 
The  right  panel  shows  the  region  allowed  for  protein  production  when  there  is  competition  for 
cellular  resources,  and  the  allowable  region  constrained  by  the  two  linear  manifolds. 
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3.  A  test  case  circuit  affected  by  competition;  A  cascade  of  transcriptional  activators 


In  the  previous  results  we  have  studied  the  case  where  two  unconnected  modules  are 
coupled  due  to  the  competition  for  the  same  resources.  Solving  that  problem  allows  us  to 
We  have  just  finished  building  a  new  synthetic  circuit  that  describes  the  case  in  which  the 
genes  competing  for  the  same  cellular  resources  are  also  transcriptionally  linked.  The 
circuit  is  composed  by  a  cascade  of  two  transcriptional  activators  that  end  up  producing  a 
fluorescent  reporter  in  the  presence  of  the  appropriate  inducers  (Fig.  4A).  The  circuit  is 
triggered  by  acyl  homoserine  lactone  (AHL)  that  binds  to  the  regulator  LuxR.  LuxR  in  turn 
activates  the  transcription  of  the  second  regulator,  NahR,  which  in  response  to  its  cognate 
inducer  salicylate  (SAL)  binds  to  the  promoter  Psal  and  promotes  the  expression  of  GFP. 


AHL  (nM) 

Fig.  4.  (A)  Schematic  representation  of  the  circuit  encoding  the  cascade  of  transcriptional 
activators.  LuxR  is  expressed  constitutively  and  produces  NahR  in  response  to  AHL.  In  the 
presence  of  SAL  NahR  triggers  the  expression  of  GFP.  (B)  Dose  response  plots  as  a  function  of 
AHL  for  different  variants  of  the  cascade  where  the  RBS  of  nahR  has  been  tuned.  The 
concentration  of  SAL  is  1  mM  in  all  the  experiments.  A  strong  RBS  for  nahR  (CAS- 1.0)  is  very 
demanding  for  ribosomes  and  synthesis  of  GFP  is  decreased  despite  producing  more  molecules  of 
the  activator.  Weaker  versions  of  nahR  RBS  (CAS-0.3  and  CAS-0.02)  display  lower  competition 
and  the  cascade  behaves  as  expected.  GFP  is  not  produced  when  the  RBS  of  nahR  is  too  weak 
(CAS-0.006)  or  there  is  no  SAL  in  the  culture  medium  (not  shown). 
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We  built  the  eascade  and  monitored  the  produetion  of  GFP  as  a  function  of  the  first  inducer 
AHL  when  the  second  was  kept  always  at  a  saturating  concentration  of  1  mM  (Fig.  4B). 
Intuitively,  an  increased  production  of  GFP  should  be  expected  when  the  concentration  of 
the  input  AHL  also  increases.  However,  the  experimental  results  show  the  opposite 
behavior  and  in  turn  less  GFP  is  formed  with  increasing  AHL.  That  observation  can  be 
explained  due  to  the  competition  of  the  two  genes,  nahR  and  GFP  for  the  same  cellular 
resources.  An  excessive  synthesis  of  NahR  uses  up  all  the  resources  and  the  production  of 
GFP  is  decreased. 

The  expected  behavior  is  nevertheless  restored  when  using  variants  of  the  circuit  with  a 
weaker  RBS  for  nahR  (CAS-0.3  and  CAS-0.06).  In  those  experiments,  the  dose  response 
plots  show  an  increased  production  of  GFP  in  response  to  increasing  concentrations  of 
AHL.  In  agreement  with  our  previous  results  using  unconnected  modules,  a  similar  output 
would  be  expected  when  reducing  the  copy  number  of  the  cascade. 

It  is  worth  noting  that  in  the  presence  of  saturating  concentrations  of  SAL  and  in  the 
absence  of  AHL  the  circuit  with  the  strongest  RBS  for  nahR  produces  a  high  amount  of 
GFP  (first  experimental  point  in  CAS- 1.0;  Fig.  4B).  That  amount  of  GFP  is  produced  as  a 
consequence  of  the  basal  transcription  of  nahR  that  may  be  originated  because  of  an 
unspecific  activation  mediated  by  LuxR  or  an  incomplete  termination  of  the  transcript 
originated  from  Plac.  In  any  case,  the  result  would  be  the  same,  a  small  amount  of  NahR 
molecules  that  can  trigger  GFP  production  in  the  absence  of  competition.  We  are  currently 
trying  to  reduce  this  basal  production  using  other  circuit  topologies  to  have  a  better 
understanding  of  the  system. 

Taken  altogether,  these  results  validate  our  framework  to  study  competition  of  different 
genes  for  the  same  cellular  resources.  In  this  particular  example,  excessive  production  of 
one  of  the  activators  has  a  consequence  the  opposite  of  the  desired  behavior. 


Current  impact: 

With  the  findings  of  this  project,  we  now  have  the  necessary  knowledge  to  determine  how 
circuits  should  be  designed  so  that  coupling  effects  due  to  unavoidable  resource  sharing  are 
minimized.  This  will  enable  the  creation  of  large  circuits  expressing  many  molecules  that 
behave  predictably  once  interacting  in  the  cell  environment. 
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